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NOMENCLATUE^ 

A CROSS SECTIONAL AREA OF THE FIN NORMAL TO xC ) 

A Surface area of the fin (m^) 

o 

Bi Biot number at the base of the fin ( h^L/k) non -<ii mens ioni 

h- Heat Transfer coefficient between the fin and the 

3 . 

surrounding fluid (V//m^.K) 

h^ Heat Transfer coefficient between the fin base 

and the fluid past the base wall (v//m JC) 

I (x) Modified Bessel function of the first kind and order n 
n' 

k Thermal conductivity of the fin material (w/m.K) 

L Length of the fin (m) 

2h^L^ 1 /2 

N A fin parameter ( ( ) 

R Radius of the spine (m) 

T Temperature (K) 

T ambient temperature (JC ) 

3 

Temperature of the fluid past the fin base wall ) 

w The width of the fin (m) 

X Dimensionless parameter (x/L) 

K^(x) Modified bessel function of the second kind and 
order n. 



'Alickness of the fin base wall 


A dimensionless parameter, corresponding to the 
thickness ratio (‘5 /l) 

Dimensionless paj^ameter (T/T^) 

Dimensionless parameter (T^/T^^) 

A dimensionless parameter defined as ( 0 -1 ) 
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SYNOPSIS 


Theoretical analyses for finding the optimal 
dimensions of longitudinal fins and spines of concave and 
convex parabolic profiles have been carried out. In the 
present analysis it has been assumed that convective boundary 
conditions exist at the fin base wall . The resistance offered 
to heat flow by the fin base wall has been accounted for. 

The analyses derive the heat transfer rate for each 

type of the fin. The optimization procedure ^then leads to a 

2h 1/2 * 

relationship between the parameters N=( ) , 6 =6/L 

hjpL ^ 

and , which must be satisfied to obtain the optimum 

dimensions for that particular fin. The relationships are 
presented in the form of graphs, Softv/are has been developed 
for design applications , Program listing is included in the 
thesis. Concave parabolic spine -is found to be the most 
efficient of the four fin types studied in this work. 



CHAPTER-I 


INTRODUCTION 


1 .1 INTRODUCTION ; 

In the design and construction of various types of ■ 
heat transfer equipment, an extended surface is attached to 
the prime surface specifically to enhance the heat transfer 
rate between a solid and an adjoining fluid. Such an extended 
surface is termed a fin. The heat is conducted through the 
material and finally dissipated to the surroundings by 
convection. An analysis of the combined conduction-convection 
systems is therefore important from the practical stand point , 

Fins offer a trouble free and compact method of enhancir 
heat transfer. Extended surface applications are numerous. 
Examples may be found in chemical, refrigeration and 
cryogenic processes ; electrical and electronic conponents , 
nuclear fuel modules, cooling fins of air cooled engines, the 
fin extensions to the tubes of radiators and other heat 
exchangers . 

A straight fin is any extended surface that is attached 
to a plane wall. It may be of uniform cross sectional area or 
it*s cross sectional area may vary with the distauice x from 
the base wall . A pin fin or spine is an extended surface of 
circular cross section. Spines can be of uniform or non 
uniform cross section. 
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Fins of different geometries and materials respond 
differently to identical and uniform heat sources or sinks . 
Important to the analyses of fin geometries are the constraints 
or the assumptions which are usually employed to define and 
limit the problem, with a view to simplifying it*s solution. 
These assumptions for conventional analysis are : 

(1 ) Steady state conduction takes place in the fin 

(2) Fin material is isotropic 

(3) Heat transfer coefficient is constant over the 
entire surface of the fin 

(^) Temperatu re of the fin surroundings is uniform 

(5) Fin thiclcness is so small compared with it*s length 
that the temperature gradients across fin thickness 
may be neglected 

(6) Temperature at the base of the fin is known and 
uniform 

(7) Resistance offered to the flow of heat by the fin 
base wall is small and may be neglected 

(8) Heat flow from the fin tip is negligible 

(9) Contact resistance between the fin and prime sxirface 
is negligible 

(10) There are no heat sources or sinks within the fin 

(11) Heat transfer to or from the fin is proportional to 
the temperature excess between the fin ar^ the 
surrounding medium. 
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For the fin to dissipate a maximum amount of heat 
for a minimum value of fin material used, it is necessary 
to use proper fins in modern designs. An optimally designed 
fin would thus minimize the cost, the space requirement 
and the weight of the system. 

One can optimize the fins in one of the following 
two ways . One is to design the shape of the fin for a maximum 
heat dissipation, when dimensions of the fin are fixed. The 
other way is to find the dimensions of the fin for maximum 
heat transfer, after the shape has been decided upon, 

1.2 LITERATURE REVIEW : 

As mentioned before, the traditional appiroach to the 
optimization of extended surfaces has been, to minimize the 
consumption of fin material for the maximisation of heat 
transmission. 

Earlier optimization procedxores carried out by Shmidt 
( 1 ) , Duff in ( 2) , Brown ( 3) , I4irray ( 4) and Gardner ( 5) were 
based on the conventional fin analysis assumptions mentioned 
in the introduction, Kern ani Kraus (6) have summarised many 
of these contributions in their book. 

The suggestion that a realistic optimization analysis 
must employ a convective boundary condition at the base of the 
fin wall and include the conduction resistance of the fin base 
wall was made by Suryanarayana (7) in 1977. This suggestion 
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amounts to the relaxation of the assumptions 6 and 7 of the 
conventional fin design analysis. Using 2-D finite difference 
method he demonstrated that the error introduced into the 
results due to these two assumptions was significant. 

An optimization analysis based on a convective 
boundary condition was carried out by Aziz (s) in 1985, But 
he did not take into account the resistance offered to the 
flow of heat by the fin base wall, 

^ An optimization analysis for various longitudinal and 
radial fins, by Relaxing the assumptions 6,7. and 8 in the 
conventional fin analysis was carried out by Alex (10) in 1986* 

1 .3 THE PRESENT WORK . 

The optimization analyses carried out so far, as 
mentioned earlier, were based on a number of assumptions . Though 
Alex (10) carried out his analysis relaxing some of those 
assumptions he left out some of the fin profiles. Not much 
work has been done to relax these assumptions at least for the 
fin profiles co^idered in this work. 

It is not realistic to assume a known base temperature, 
since in an actual situation of heat dissipation by a fin 
the heat transfer is governed by cony<ection. For a more accurate 
analysis one must apply a convective boundasr/ condition on 
the plane side of the fin base wall. 
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The conventional analysts of the fins also ignores 
the resistance offered by the fin base wall. This is also 
not realistic. So in this analysis we have taken into account 
the resistance offered by the fin base wall to the flow of 
heat. From the graphs for the optimum dimensions of fins 
( Fig .Nos .5, 7,9,11 )» we see that the inclusion of the 

conduction resistance offered by the fin base wall, in the 

-1 

analysis, lowers the value of Bi and N for constant values of 

-1 

N and respectively. 

This thesis includes the results of the optimization 
analyses of longitudinal fins of convex and concave parabolic 
profiles and spines of concave and convex parabolic profiles, wi‘ 
the assumption numbers 6 and 7 of the conventional analysis 
relax:ed . 

For each type of the fin the basic differential 
equations governing the heat transfer -in the fin and fin base 
wall are solved simultaneously, to yield an expression for 
temperature excess. Next the correct temperatxire distribution is 
obtained using the appropriate boundary conditions. The 
equation for the heat dissipation by the fin is then obtained 
from the Fourier's irate equation. 

The optimization procedure employed is to express 
the dissipation, q, as a function of, w, the width of the fin 
in case of longitudinal fins of concave and convex parabolic 
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profiles, and then to apply the condition ^ = 0, for 
maximizing the heat transfer. In case of spines of concave 
and convex, parabolic profiles, q, is expressed as a function 
of R, the radius of the spine and then the condition ^ = 0 
is applied for maximizing the heat transfer. 

For each type of fin or spine, the above optimization 

procedure leads to a relationship between the fin parameter, N, 

* 

and the biot number, Bi , and the thickness ratio, d , that 
must be satisfied for optimum conditions. The variation of, 

Bi, with , N, under different conditions are then plotted. 

Since we have not assumed a known base temperature an 
iterative solution is required for the convective boundary 
condition. 

For a desired heat dissipation rate, the fin design 
procedure would be as follows ; delect a trial value of, L, 
the length of the fin, which together with other data fixes, 

Bi, andjd* . Using the Relationship between N, Bi and 6- , N 
value can be found out. Now we calculate the heat dissipation, 
q, using the appropriate expression for the desired fin type 
and compare it with the desired heat dissipation. We repeat the 
procedure by incrementing the length till the values match. 

From the value of N, at which desired, q, is achieved we find ou 
the radius or width of the fin under consideration. 
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CHa RL’ER- II 


I'l/i'HSM.vTICAL POi^''IUL^\TIOK AIJD OPTIMIZATION ANALYSIS 


2.1 LONGITUDINAL PIN OF CONCAVE PAR:iBOLIC PROFILE 


Derivation of the governing differential equation for the 
flow of heat through the fin : 


The general differential equation for any longitudinal 


fin is 


,2m dfo(x) 2h_ 

^ 2 '^) ^ * 2 -Ir- i - 0 


... ( 2 . 1.1 


The profile function f 2 (x) in case of the concave parabolic 
fin shown in Fig. 1, is 


fjCx) = § (g)2 


( 2 . 1,2 


wx 


dx 


... (2.1.5 


Substituting the Eqs, (2.1,2) and (2,1,3) in Eq. (2.1.1), 
we get 


^ * 2-i- ^ ° ....(2.1.4 


dx‘ 


Eq. (2.1.4) is then the governing differential Equation for 


the concave parabolic fin 
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The differential equation governing the flow of heat through 
the fin base wall is 

d^T/dx^ =0 *••• (2.1.5) 

BOUNDARY CONDITIONS : 

riT 

1. At the tip of the fin (x=0), = 0 

2. At the wall surface where x = 6 +l, 

= hf (Tf-T) 

3. At the base of the fin i.e, at x = L, 

dT. dT^ 
dx ~ dx 


Nondimens ionalizing the governing equations for the fin, the 
base wall, and the boundary conditions by using the following 
expressions : 

T„ _p 2h 

0= T/Tg^, X= ^ ^ , N^= kw 

dT '^a d© d^ _ d^© 

a5; = irHz = ^7^ 


(1) For the fin, we get 

^ + 2X^ - N^(© -1) = 0 

dQ 

At X=0 i.e. at the tip of the fin ^ = 

(2) For the fin base wall 


0 


d^& 

dX^ 


=0 
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At X= 1+6 i.e. at the surface of the wall 

^ ©) 

(3) At X= 1 i.e, at the base of the fin 

d©. d©p 

®1 " ®2 37 “ = 37 “ 

To solve the above equations let us assume that 


©- 1 = 1 ' 


Then the equations and the boundary conditions become 


(1) For the fin 




.... ( 2 . 1 . 6 ) 


"the "tip of "tho fin i #6 • sft X=0 f = 0 


(2.1.7) 


(2) For the fin base wall 
_ 0 


( 2 , 1 . 8 ) 


At the surface of the wall i.e, at X = 


37 = B^(©^-1-1') 


( 2 . 1 . 9 ) 


(3) At the base of the fin i.e. at X =1 


^1 = ' 5^2 


.... ( 2 . 1 . 10 ) 


d^^ dl*^ 

3r‘ “ 3r" 


( 2 . 1 . 11 ) 
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The boundary conditions (2,1.10) and (2.1,11) are common to 
both the Eqs. (2.1.6) and (2.1.8). The solution of Eq. (2.1.6) 
is obtained as follows . 

Let X= e^ or v=lnX then 

dW _ d'lf dv _ 1 dW 
~ Y. 3v 




Now Equation (2.1 .6) can be written as 

....( 2 . 1 . 12 ) 

The solution of Eq. (2.1,12) is 

P.v P^v 

'1' = C.J e +02 0 

In terms of independent variable X 

P. P? 

T= c., X ' + C 2 X .... (2.1.13) 

1 /2 

where P., = -1/2 + 1/2 (1+4N^) 

1 /2 

= -1/2 - 1/2 (1+4N^) 

It can be observed that at X=0, the temperature excess will 
be un-bounded unless C2=0 in Eq, (2.1.13), 

T = c.^ X 


• • « • 


(2.1.14) 
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The solution of Eq. (2.1.8) is 

'i’ = AX +B (2.1.15) 

Using the boundary condition (2.1.9) in Eq, (2.1.15), 
v/e get 

A = B^(©^-1-A(1+5*)_b) 
or B = ( ©^—1 ) — 1 \ ( 1 /B j^+1 +6 ) 

'i'=AX+(©^-1) -A(1 /B^+1+6*) (2.1.16) 

Similarly using the boundary condition (2.1,10) in Eqs, 

(2,1,14) and (2.1.16), we obtain 

c^= (©^-1) -A(l/B^+6*) (2.1.17) 

Substituting the boundary condition (2,1,11) in Eqs. (2.1,14) 
and (2.1.16), we have 

A = P^c-j .... (2.1 .18) 

Also substituting Eq,(2.1 ,18) in Eq, (2,1.17) yields 
c^ = (©^-1) - P^c^(l/B^+6*) 

( ©^ “1 ) 

c-= 

' 1+P^(l/B^+6 ) 

Substituting the value of c^ in Eq, (2.1,14) gives 

P'1 * 

$=(©^-1)X 7(1+Pi(1/B^+<S )) .... (2.1.19) 

P^X^I (©^-1 ) 

1+P^(l/Bj^+ <5 ) 


. d'F 
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1^P^(l/B.+6*) 


The heat transferred through the base of the fin is. 


q = KA (^) 


X=L 


q = 


KW (©f-1)P^ 


1+P^(l/B.+6 ) 


( 2 . 1 . 20 ) 


To optimize the heat transfer let us assume that liie volume of 
the fin, 1/3 V/L=V (a constant). 


‘The optimization procedure is to express, q, as a 

function of, w, alone and then applying the condition ^ = 0. 

The application of this condition yields the following 

.* 

relationship between B^,o and N, 

1 /2 

1 /R = ((1+4N^) ~N^-1) _5* 

i (1+ An ^)'' /^-1 )^ 

(For derivation please refer to appendix A) 


2.2 LONGITUDINAL FIN OF CONVEX PARilBOLIC PROFILE : 

Derivation of the governing differential equation : 

The general differential equation for any longitudinal 


fin is 


dfp(x) 

2f2(x) ^ +2 “■Jx"* 35E 

ax 


d^T 


2^a 

IT 


(T-T^)= 0 


• • • ♦ 


( 2 . 2 . 1 ) 
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The profile function f 2 (x) in case of the convex parabolic 
fin, shown in Fig, 2, is 


= I (f ) 


1/2 


( 2 . 2 . 2 ) 


df^Cx) 

3x 


w 


4 ( Lx) 


1/2 


(2.2.3) 


Substituting (2.2,2) and (2.2.3) in equation (2.2.1)^ we get 

( I) (T-T^) =0 (2.2.4) 

1/2 

where m = (2 ha/kw.) 

Equation (2.2,4) is the governing differential 
equation for the convex parabolic fin. 

The differential equation governing the flow of heat 
through the fin base wall is , 


d^T _ Q 

dx^ 


(2.2.5) 


dT 


BOUNDARY CONDITIONS : 

1 . At the tip of the fin When, x * ^ 

2. At the wall surface when x = L+6 , 

^ ix“^f 

At the base of the fin i ,e . at x = L , 
dT^ dT^ 

T^=T 2 and 33 ^ = ^ 


3 . 



Nondimensionalizing the governing equations for the fin, 
the base v/all, and the boundary conditions by using the 
following expressions . 

& = X= 0^= ^ 
dT "^a d© d^T "^a d^© 

Then, 

(1 ) For the fin , we get 

1/2 2 /2 

d^© ^ X 1 d© ,^2/0 IN n 

^ ^JT72 33? 

d© 

At X=0 i.e, at the tip of the fin ^ = 0 

(2) For the wall 



A X= 1+6*, i.e. at the surface of the wall 

if = ©) 

j 

At X=1 i.e, at the base of the fin 

d©^ d©2 

©1= ®2 33r = 33r 


To solve the above equations let us assume that © -1 =?? 
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Then the equations and the boundary conditions become 
(1 ) For the fin 




I , d.'i^ TvT^m 


At the tip of the fin i.e, at X=0 ^ 


= 0 


(2) For the wall 


.... ( 2 . 2 . 6 ) 
.... (2.2.7) 


dh!_ 


.... ( 2 . 2 . 8 ) 


At the surface of the wall i.e, at X= 1+6 


II = B^(©jp-1-¥) 


.... (2.2.9) 


At the base of the fin i.e. at X=1 

^^=^2 *•" ( 2 . 2 . 10 ) 

df. d?5 

" 3r" •••• (2.2.11) 


The boundary conditions (2.2.10) and (2.2.11) are common to 
both the Eqs, (2.2.6) and (2.2.8). 

The solution of Eq. (2.2.6) is. 




1 /4 3/4 

=X C^I^/3(4/3N X )+ 


Col 


2-^-1 /3 


(4/3N 




• • ♦ • 


(2.2.12) 



The use of boundary condition (2,2.7) in equation 
(2.2,12) requires multiplication of each of the terras of the 
infinite series expansion of I.jy5(4/3N X^^^) and X^^ ) 

by followed by a terra-by-tenn differentiation. When this 

is performed the temi involving (^) (x"^ ^^(4/3N X^^^) ) 

becomes unbounded at X = 0 which necessitates that C.j=0, 

Therefore T= x"* /5(4/3N X^'^'^) .... (2,2.13) 

The solution of equation (4,2,8) is 

'!'= AX +B .... (2.2.14) 

Using the boundary condition (2,2.9) in (2,2,14), we obtain 
A=B^(©f-1-A(1+6*)-B) 
or A(1 /b^+1+6*) = (9^-1 )_B 

T=AX+(6^-1)-A(1/B^+1+6*) .... (2.2.15)’ 

i 

Using the boundary condition (2.2,10) in Eqs. (2,2wl3) and 
(2,2.15), we have 

C2l_i/3(4/3N) = (©f-1)-A(l/B^+6*) .... (2.2,16) 

Similarly substituting the Boundary Condition (2.2.11) in 

Eqs. (2.2.13) and (2.2.15) yields 

A= Nl2/3(4/3N) .... (2,2.17) 

Also by substituting (2,2.17) in (2.2.16), we get 
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or 


^2^-1 N)= ( 9 ^- 1 ) -02^2/3(4/3 N) (l/B^+6*) 

C2(I_^ 73(4/3 N)+N 12/3(4/3 N) (l/Bj,+6*)) =©^-1 

&.-1 

Q - ^ 

^ 1 - 1 73(4/3 N)+Nl 2 / 3 ( 4/3 N) (l/B^+6 ) 


Substituting the value of C2 in Eq. (2.2.13), gives us 


x''/^( 0^-1)I_^/3(4/3 N X^/^) 
I_1/3(4/3 N)+Nl2/3(4/3 N) (I/B^+6''') 


(2.2.18) 


dW 


(Q^-1) (Nl2/3(4/3 N X^^^) ) 

1.1 73 ( 4/3 N)+Nl2/3(4/3 N) (1/B^+6*) 


(2.2.19) 


The heat transferred through the base of the fin is 

T, 

C 

L 


vA / dT \ T^A a 

q = KA (37)^^£ KA ^ 


kwT, 


or 


q = 


(9^-1)(Nl2/3(4/3 N)) 


73(4/3 N)+Nl2/3(4/3 N)(l/Bj_+6 ) 

.... ( 2 . 2 . 20 ) 


To optimize the heat transfer let us assume that the 
volume of the fin | wL=V (a constant). 
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The optimization procedure is to express, q, as a 
function of, w, alone and then applying the condition, ^ = 0, 

The application of this condition yields the following 
relationship between B^, 6 and N, 

^ ^ 2N^(I_^/3(4/3 N))^-(l2/3(4/3 N) )^-Nl2/3(4/3 N)I_^ ^3(4/3 N) 

^ ^ N^(l2/3(4/3 N ) 

( For derivation please refer to Appendix B) 


2.3 CONCAVE PARABOLIC SPINE 


Derivation of the governing differential equation: 

The generalized differential equation for any spine is given by 


2 .2, 


(fjCx)) ^ ^ (f2(x) -IT ^2(^5 (T-T^)=0 

.... (2.3.1) 

The profile function ^2^^^ concave parabolic 

spine, shown in Fig. 3, is 

- ....(2.3.2) 


>2 dT ^a 


f2(^) “ 


h. ^2^^^ 


D 


(2.3.3) 


Substituting Eqs . (2.3.2) and (2.3.3) in Eq. (2.3.1), we get 

. 4x3 i - ^ x2(T-T,)= 0 

dx 


.... (2.3.4) 
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Eq. (2,3.4) is the governing differential equation for a 
concave parabolic spine . 

The differential equation governing the flow of heat in the 
fin base wall is 



.... (2.3.5) 


boundary conditions : 

1 . At the tip of the fin where x = 0, ^ = 0 

2. At the wall inner surface, where x =L+6, 

‘'ii = hf(Tf-T) 

3. At the base of the fin i.e. at x=L, 


dT. dTp 
T^=T2 and ^ ^ 


Nondimensionalizing the governing equations and boundary 
conditions : 

^ T-f o 2h^L^ 

Let 0=T/Tg^, X= ^ ,0^= Nr= 

3 . 

Then 

(1 ) For the fin, we get 

^ + 4 X ^ ^ - N ^ X ^(0 - 1)=0 
dX^ 

d© 

At X=0 i.e. at the tip of the fin = 0 

(2) For the fin base wall 



0 
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At X= 1+6 ^ i.e, at the surface of the wall 

II 

At X=1 , I.e, at the base of the fin 

d&. d&p 

©2 and ^ ^ 

To solve the above equations let us assume 
T =9 -1 


Then the equations and boundary conditions become 

(1 ) For the fin 


^ +4x^ Iy -N^X^ = 0 


d? 


, . , (2,3,6) 


At the tip of the fin where X=0 
d'F 


d^ 


=0 


(2.3,7) 


(2) For the fin base wall 



.... (2.3.8) 


At X= 



i ,e , at the inner surface of the base wall 


g = B.CSf-l-f) 

(3) At the base of idie fin, where X=1 


d?^ d¥2 

37“ “ 


(2.3.9) 


(2.3.10) 

.... (2.3.11) 
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The solution of Eq. (2,3.6) is 

(2.3.12) 

1 /2 

where =(-^2+1/2 (9+4N^) ) 

1 /2 

=(-3/2-l/2 (9+4N^) ) 

The value of, C 2 , in Eq. (2.3,12) must equal zero for obtaining 
a finite temperature excess at the tip of the fin (where X=0) 


'I' = C^X ' 

The solution of Eq, (2,3,8) is 
=AX+B 


.... (2.3.13) 


.... (2.3.14) 


Substituting the boundary condition (2.3.9) in Eq. (2.3,14), 
yields 

B= (©£-1)-A(1/B^+6*+1) 

/, 'i'=AX+(©^-1 )-A(1 /B^+4 +1) •••• (2.3.15) 

Applying the boundary condition (2.3.10) to Eqs. (2.3,13) 

« 

and (2.3.15) , we get 

C^= (9^-l)-A(l/B^+6*) 

Applying the boundary condition (2.3.11) to Eqs. (2.3.13) 
and (2.3.15), we obtain 


.... (2.3.16) 


A= P^ 


.... (2.3.17) 
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oubstituting the value of Eq, (2.3.17) in the Eq. (2.3.16), 


gives us 

C = ^ 

^ (1+P^(1/B.+6*) ) 


Substituting the value of, , in equation (2.3.13), yields 


(9„-1)X 

^ 

(1+P^(l/B^+6 ) ) 
P.-1 

^ (1+P^(l/Bj.+6'^) ) 


(2.3.18) 


.... (2.3.19) 


Theheat transferred through the base of the fin is. 


q 



KnR^ P-,(©f-'l) 

^ ’(1+P^(1/B^+6*)7 

.... (2.3.20) 

The procedure for optimization involves expressing q 

as a function of P alone and applying the condition =0- 

The optimization analysis yields the following relationship 
* 

between, ^ and N* 
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1/2 1 /2 

^ 5N^+2(9+4n^) (3(9+^N^) -1) 

^ 2(9+An^)^'^2(^(9+4N^)'’^^-3/2)^ 

( For derivation please refer to Appendix-C) 


2.4 CONVEX PARABOLIC SPINE : 


Derivation of the governing differential equation: 
The profile function, of a convex parabolic spine 
(shown in Fig. 4) is given by 





(2.4.1) 


df2(x) ^ ^1/2 


.... (2.4.2) 


VVe have the general differential equation for a spine given by, 

(fjCx) ^ - I; (f2(x) i - ^ fat’') f-0 t2.^.3) 


Substituting Eqs. (2.4.1) and (2.4.2) in Eq. (2.4.3), 
we get 

1/2 

d^T ' dT (T T ^ 

^7-7 ^ 35E — m — (T-T^)- 0 

dx 


(2.4.4)' 


Equation (2.4.4) is the governing differential equation for the 
flow of heat in a convex parabolic spine. 

The equation governing the flow of heat in the 


spine base wall is, 
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BOUND.- ;RY CONDITIONS : 


1 . 
2 . 

3 . 


At the tip of the spine (x=0), ^ =0 

At the wall inner surface (x=L+6), 
k dT/dx=h^(T^-T) 


At the base of the spine (x=L), 

dT. dT^ 

^1=T2 H5r = dT" 


.... ( 2 . 4 . 5 ) 


Nondimens ionalization of goveminf^ equations and boundary 
conditions : 


Let » =T/Tg^, 



,©^= wr“» 

a 



2h 

a 

1^ 


Substituting the above values in the differential equations 
and boundary conditions, we get 


1) 


At 


For the spine 




dX 


. d& 


X=0, If - 0 


1 /2 

-N^(X) (O -1)=0 


2) For the spine base wall : 


d2g 


= 0 


dX‘ 

-* d© 


At X=1+<5 , ^ ©) 
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3) At X-1 , i ,e . at the spine base 
®1"®2 dl" = 35r 

Let 'S= © -1 

Then, the above equations and boundary conditions become 


(1 ) For the spine 


-if * s =0 


(2.4.6) 


At X=0 i.e, at the spine tip 



(2) For the spine base wall 



.... (2.4.7) 


.... (2.4.8) 


At X= 1+6*, 11= .... (2.4.9) 

(3) At X=1 

...» (2.4,10) 

d^. d'Pp 

and .••• (2.4,11) 

The boundary conditions (2,4,10) and (2,4,11) are 
common to both the equations (2.4,6) and (2.4,8) , The solution 
of equation (2.4.6) is, 

'}'=C^(4/3 N X^/^)+C2k^(4/3 N X^^^) 


(2.4.12) 
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For satisfying the boundary condition (2.4,7) i-.e, 
to maintain a finite temperature excess at X= 0 ,C 2 must be 
zero since is unbounded, 

3/4 

^=C^I^(4/3 NX) (2,4,13) 

'Ilie solution of Eq, (2,4,8) is 

'I'=AX+B (2.4,14) 

Using the boundary condition (2,4.9), Sq. (2.4.14) can be 
v/ritten as 

B= (9^-1) -A (1/B^+1+6*) 

Substituting the value of B in Eq. (2.4,14), we get 

W=AX+(9^-1 )-A(1/B^+1+6*) (2.4.15) 

n.pplying the boundary condition (4,2,10) to both the Eqs, 
(2.4.13) and (2.4.15), we obtain 

Cl Iq(4/3 N)=(©^-1)-A(1/B^+6'‘) (2.4.16) 

Using the boundary condition (2.4,11) in the Eqs, (2,4.13) and 
(2,4.15) yields, 

A=C^N.I^(4/3 N) (2.4,17) 

Substituting Eq. (2.4.17) in Eq, (2.4.16), 
we can find out the value of 
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C.= 


C^.I^(4/3 N) =(9^.1)_C^N. I^(4/3 N). k 
(©f-l) 

IqC 4/3 N}+NkI^(4/3 N) 


i -H- 

v/here k= ^ +6 
i 


oubs'ti'tuting the value of in (2,4,13), gives us 

3/4 

(4/3 N X ) 


'?= 


T7(^^73Tr)Tm:;W3r¥r 


-1/4 3/4 

X NI^(4/3NX )(9^-1) 

33r Iq( 4/3 N)+N,k.I^(4/3 X') 

The heat transferred through the spine base is 


q = 


X=L 


K;\T 


a 


X=1 


(2.4.18) 


(2,4.19) 


a = 


knR^T^(©j^-l)N.I^(4/3 N) 
l(i:"(Wn)4i, (4/5-TJVK.i^. )’ 


(2.4.20) 


Applying the optimization criterion, ^ =0 leads to 

* 

the following relationship betv/een B^,<3 and N 

6 5 2lo(V3N) 

(I,(4/3n))^ -5-TII,(V3 M) 

(For derivation please refer to appendix-D) 
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CHAPTER-III 

RESULTS AND DISCUSSIQI^IS 
3.1 RESULTS A^3D DISCUSSIONS : 

3.11 . LONGITUDINAL FINS OF CONCAVE PAR.-.BOLIC PROFILE : 

Figure (5) depicts the optimum dimensions for concave 
parabolic fins. It represents the relationship, between N and 
Bi that must be satisfied by a fin to transfer the maximum 
amount of heat . It is seen from the figure that as N 

* 

increases the value of 1 /Bi decreases, for a particular o 
value. For a higher value of N, a higher Bi is necessary to "bransl 
the same amount of heat. Therefore 1 /Bi decreases as N increases, 
as b* increases the curves become more steeper and for a constan't 
value of N, 1 /Bi decreases . v^hen the thicknes of the fin base 
wall is more, the resistance offered to the flow of heat by 
the wall is greater and so heat transferred is less. To 
compensate this reduction in heat transmission a higher Bi 
value is required . 

Figure (6) shows the variation of the fin efficiency (h) 
with the parameter N for constant values of (1/Bi+6 ). As the 
value of N increases the value of the efficiency decreases , 

It is observed from the figure that as (1/Bi+6 ) increases, 
the curves become steeper as they shift diagonally to the left 
and for a constant value of N, the efficiency decreases. 
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'H' 

The top most curve, in v/iiich the value of (l/B^+6 ) is 
equal to zero, represents the simplest case when the base 
temperature of the fin and the fluid temperature are equal. 

This corresponds to the situation v;here the heat transfer 
coefficient, h^, between the fluid and the fin base wall is 
infinite, and the fin base v/all resistance is zero. For higher 
values of (I/B^+S ), h^. becomes finite and therefore introduces a 
resistance to the flov; of heat from the fluid to the fin base 
wall, in addition to the resistance of fin base wall itself, 
v/hich results in a decreased efficiency. 

3.1.2 LONGITUDINy'iL FIN OF CQWEX R.th.BOLIC PROFILE. 

Figure (7) is the plot of the optimum dimensions for 
the convex parabolic fins. It is a graph of 1 /B^ vs N for various 
values of 6 . As N increases 1 /B^ decreases. For higher values 
of 0 , the decrease in 1 /B^ becomes steeper and shifts 
towards lower N. 

Figure (8) gives the variation of the efficiency of 

convex parabolic fins with M. It is seen that as N increases 

x- * , 

the efficiency decreases for various a values. As 6 increases 
the curves become steeper and shift diagonally towards lower 
values of N and for a constant N the efficiency suffers a 
reduction in magnitude. This is, as before, due to the 
reduction in heat transfer ^caused by the increased resistance 
of the fin base Jail. • ■ 
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3 .1 ^ CONCAVE PARABOLIC SPINES : 

Figure (9) shows the optiraum dimensions for concave 

parabolic spines, for various 6* values. This 1 /Bi vs N plot 

represents the relationship between 1 /Bi and N that has to-be 

satisfied to transfer the maximum amount of heat through the spine. 

As N increases 1 /Bi decreases, for a particular 6 value. Curves 

* 

become steeper as o increases giving rise to a larger drop 
in 1 /Bi with N . 

Figure (10) presents the change in the efficiency of 
concave parabolic spines with N, for different values of 
(1/Bi +‘^ ). Efficiency is seen as usual to decrease with 
increasing N, for a particular (1/Bi+6 ) value. It can also be 
observed from the figure that a decrease in efficiency occurs 
as a result of increasing (1/Bi+6 ) value, for a constant 
value of N. 

3^.4 CONVEX PARABOLIC SPINES ; 

Figure (11) exhibits the effect of 6 on the optimum 

dimensions of a convex parabolic spine. It is a graph of 

* ♦ 

1 /Bi Vs . N, for different 6 values . For a given value of o , 

1 /Bi is seen to decrease with increasing N. An increase in the 

value of 6* marks a decrease in 1 /Bi for a constant value 


of N. 
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Figure (12) is the plot of efficiency Vs. N, As is 
the case before, efficiency decreases with increasing N, for a 
constant value of (l/Bj^+6*), As (1/^^+6 ) increases, the 
efficiency decreases for a particular value of N, 

Figures 13,14,15 and 16 give the optimum dimensions 
of various types of fins for 6 =0, 1,2, and 3^ respectively , 

These are plots of 1 /B^ Vs. N. For a constant value of N less 
than, 1 , the concave parabolic spine gives a higher value of 
1/B^, compared to all other analysed fin types. This means 
that for a given N value below 1 , the optimal length required 
for the concave parabolic spine is going to be the least 
compared to the other analysed fin configurations which yield 
lov;er 1 /B^ values . 

Figures 17,18,19 and 20 depict the variation of the 
efficiency of various fins with the fin parameter N, for 
(1/B^+6*)= 0,1,2, and 5, respectively. It is seen that the 
concave parabolic spine has the highest efficiency and concave 
parabolic fin the lowest. The efficiency of other fins lies 
in between these two. The efficiency of a particular fin 
geometry is seen to increase with decreasing N values. 

Figure 21 gives the optimum dimensions of various fins 
analysed In the present work and in a previous one (Ref .10). 

It is a plot of 1/Bj^ Vs N for 6*=0. 
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, * 

Figure 22 is a plot of efficiency Vs. N for 1 =0 
for various fin profiles, present and previous (Ref .10). 

As can be seen from the graph concave parabolic spine is -the 
most efficient of all the fin types considered. 
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CONCLUSIONS 


main conclusions DFL\WN from this study are as follows : 


1 , The efficiency of all the fins investigated herein 

decreases as the fin parameter N increases. 

2. For a constant value of fin parameter N, the 

* 

efficiency decreases as the value of (1 /Bi+<5 ) 
increases . 

3 . The value of 1 /Bi decreases as the fin parameter, N, 
increases , 

4. For a constant value of N, 1/Bi decreases as 6^ increases 

5 . Out of the four types of fins analysed the concave 

is 

parabolic spine/the most efficient and the concave 
parabolic fin the least efficient. 

6 . Under the optimum condition , the concave parabolic 
spine yields a maximum value of 1/Bi, for a constant 
value of N, while the convex parabolic spine gives a 
minimum value. 

parainG*t6rs Bi and 6 hatvo a. considerabla effsct on 
*bhs diniBns ions oi "tbG fins • 
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3.3 SUGGESTIONS FOR FUTURE V/QRlx . 

A similar analysis can be carried out for other 
types of fins like, longitudinal fins of trapezoidal profile, 
truncated conical spine, radial fins of concave parabolic, 
convex parabolic, hyperbolic, triangular profiles etc. 

In the present study constant fluid and fin material 
properties have been assumed. Thermal dependence of these 
fluid and fin material properties could also be considered 
in a future optimization analysis. 
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APPENDIX -A 

From Eq. (2.1.20) on page (12) 
V/e see that 

P.w^ 

q a ar- 

1+P^(l/B.+6 ) 


The optimization criterion is |^ = 0. 

— i /P 

(1+P^ (1 /Bj^+6 ) ) .(P^ .2w+w^((4N^+1 ) 2N(^“) )) 

P 1/B.+6* „ ^ _'l/2 

= P^w^(P^( i )+(1/B^+6 )(4N^+1) 2N(^) ) 

(2.1.21) 


2 iShV^ 


where V= -^Lw 


. dN -3N 

aw = 


* 



d f kw 6w 

aw'ih^V 37 

_ k . 1 ,5* 

- T3CT 37^" F3? w" 

=■(—4 — ) 

LHS. of Eq. (2.1.21) = 

—1 /2 -2. 

(1-fP^ (1 /Bj^+6 *) ) (P^2w+w^(1+4N^) C-^^) ) 

= (1+P^(1/B^+6*))(P^2w+ ^(-3n2(1+4n2)-W2)) 
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f 2 

= f1+P^ (1 /B^+6 ))(w(2P^-3N^(1+4n^) )) 

RHS of Eq . (2.1 .21 ) 

2 /p 

Pl^ (P'l(l/Bi+6*) + (l/BjL+6*).(1+4N^) (-3N^)) 


•H" "“*1 

= P ^ w ( l / B ^+ 6 *)( p ^> 3 N ^( 1 + 4 N 2 y j 


Eq. (2.1,21) can now be written as, 

(1+P^(1/B^+6 ) )w (2P^ -3N^(1+4n^) ) 

—1 /2 

= P^ w(1/B^+6*)(p^-3N^(1+4N^) ) 

•“1 /2 

.*. 1+P^ (1 /B.+6'')(2P^_3n2(1+4n2) =p^ (1/B^+6*) 

-1 /2 

x(P^-3N^(1+4n^) ) 

^ /2 -*‘1 j ' 2 m 

2P^+2p2(1 /B.+6*)-3N^(1+4N^) -P^ (1 /B^+6*) 31^(1+41^) 

—1 /2 

= P^(1 /B^+6*) _p^ (1 /B^+6*)3N^(1+4 n^) 


—1 /2 

.*. 2F^+P^(1/B^+6*)-3N^(1+4N^) =0 

-1 /2 

.•. P^(1/^^+6*)=3N^(1+4N^) -2P^ 


(1/B.+6*) = 


-1 /2 

3N^(1+4N^) -2P. 


.1/2 

Replacing P^ by the expression, P-^= -1 /2+1 /2 (1+4N^) in 
the above equation yields 
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(-I/2+I/2 (1+4n^)'^'^ 


1 /2 

5N^+ (1+4n^) - (1+4N^) 

T T72 — ^ Z 

(1+4N^)/^(1 //j+1/4 (4N^+) -1/2(1+4N^) ) 


4 (3N^-4n^+1+(1+4N^) 

' i/d — : 

(1+4n ) ((1+4N‘^) -1) 


•*• h: = 


4( (1+4n^)^^-N^-1 ) 

^ ; t 72 — : 

(1+4N^) ((1+4n^) -1) 


4( (1+4n^) 

— '175 •■-' rr/rn 

(1+4N^) ( (1+4N^) -1) 
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Ai^Pi'HDIK-B 

From Eq. (2,2,20) on page number (17) we see that 

w^(NI^ /,(4/3 N) 

q a ^ ^ ^ 

I_^/3(4/3 N)+Hl2/5(^^/3 ^)(l/H^+b ) 

Application of the optimization condition -0 yields, 

(I_^/2(4/3 N)+Nl2/3(4/3 N) (1 /B.+6*) )|_(w2nI2/3C4/3 N)) 

= n))^(I^^^^^^/3 n)+NI2/3(4/3 N)(1/B^+6*)) 

( 2 . 2 . 21 ) 

where 

‘ ■’^-1/3^'*^^ ^^2/3^'''^^ 

= |!! 1/2Nl2/3('./3 N)-I_^y3{4/3 N) 
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I; I.i/3(V3n) =I_^/3(V3 N)+ 1 I , (4/3 N) 

33N 

3lU 

- - Ij,/3(V3 N)45,i.^/ 3(4/3 N)ifi 

" -l2/3(V3 N)- V-1/3 

(4/3 N) |5 

° 5H ^-1/3<''''2 N)-l2/3('t/3 N)P 

3w ^-1/3^^^^ '2w ^-1/3^^^^ 

substituting the above values into Eq. (2.2.1) we get 
L.H.S. of (2.2.21) 

^^-1/3^^/^ N)+r^^/3(4/3 N)(l/B.+6*))(Ni2/3(4/3 N) 2w+w^ 

= I_^/3(4/3 N)+Nl2/5(4/3 N)(1 /B^+6*)2 w Nl2/3(4/3 N)- 

l2^3(4/3N)+Nwl2^3(4/3 N)-2N^ ^ 3 ( 4/3 N) 

= I_1/3(4/3 N)+Nl2/3(4/3N)(l/B^+6*)_3/2Nw 12/3(4/3 N)- 

21^1.1/3(4/3 N) 


41 


= n)i^/3(4/3 ^^2 

(^-1/3(^/3 M)) -2KfA, 12/3(4/3 N)I_^/3(4/3 n) 

41 S0 R.H.S. of Eq. ( 2 . 2 , 21 ) 

( (Nl2/3(4/3 N))(^(i_.|^^( 4/3 N)) + (1/B3t6*)(i^^^(4/5 N)(J^) 

u. 2N^/1 ^ . 


w N)) 

- 1*12/3(4/3 N) 

= w2(®2/3(V3 N))*^i_^^^(4/3 N). |Ni^^3(4/3 N)-(1/b^.6*) 

I 2 / 3 (V 3 N)i . ilL£Lx^/^(4/3 N)-2 n2 '2^> 

Ii- 3(*‘/3 N)+ — ^m 3/3(4/3 N) 

= l 2 / 3 (V 3 „,) (1/2 11^ ^^^^,4/3 

1 , 2N^(1/B.+6*) 

3 ? I-i/3(''/3 n) i I..|/3(4/3 N)) 

= 1/2(1 /B3+6*)n2w(I2/3(4/3 n))^ 

* N)I.^/3(4/3 n)-2n\ 12 / 3 ( 4/3 N) l.i/3(4/3N)(1/i,^.S*) 

iiiquating L.H.S. and R.H.S. of Eq, (2.2.21) yields 
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3/2rKv 11 ) 12 / 3 (^ 1/3 N)+ |(l/B^+6*)N^iv(l2/3(4/3 N) 

-2N^(I_^/3(4/3 N))^-2(1/b^+6'')n\ l2/3(V3 N)I_^/3(4/3 N) 

^ ^2i ^ 

= 1/2(1 /Bj_+6 )N^w(l2/3(4/3 I'])) -211^(12/3(4/3 M)) 

+ l2/3^^/^ ^'01.1/3(4/3 M)-(2n\(l/B.+6*)i2/3(4/3 N) 

+I_l/ 3 (^/^ ‘0 ) ( 2 . 2 . 22 ) 

Rearranging the terms In Eq. (2.2.22), gives 

fiv; 12/3^^^^^ N)+( 1 /Bj^+ 6 *)N^(l^^ 3 ( 4/3 n) ) 

2 ^ 2 

- 2 n 2 w CCi.i/ 3 ( 4/3 n)) -( 12 / 3 ( 4/3 ; 0 ) ) = 0 

2 

(l/B.+6*)N^(l2/3C4/3 N) ) 


= 2n2w(I_^/3(4/3 N)) .( 12 / 3 ( 4/3 I!) -Nwl2/3(4/3 N)I_^/3(4/3 ^J) 




h6* N)) '-( 12 / 3 ( 4/3 N)) -Nl2/3(4/3 N)I.^ /3(4/3N ) 

N^(l2/3(^/^ N) ) 


1 

17 


2N' 


'((I,1/3(4/3 r0)^-(l2/3(4/3 iQ) >-Nl2/3(4/3 N)l^/3(4gN) ^ 
N^(l2/3(^/3l0 ) 



43 


APPENDIK-C 


From equation (2.3.20) on page number (22) 


it 


qa 


( 1 + P ^( l / B ^+6 ) 

- applying the condition ^ =0, we get 
( 1 + P ^ k )^( R ^ P ^ )= R ^ P ^ ( 1 + p ^ k ) 

v/here P^ = (4* ^(9+4N2)''''^)'k= (^6*) 

L.H.O. of Eq. (2.3.21) 

,4 d 


1/2 


= (1 +P^k)( 4R^P^+H^ .3/2+1/2(9+41^) )) 

r 2 ‘ —I /2 

= (1+P^k)(4R^P^+R^(- ^(9+41^)” )) 

Since 


= (1 + P^k)(4R^P^-5NV(9+4N^) ' ) 

—1 /2 —1 /2 

= 4R^P^+4R^P^k-5N^R^(9+4N^) -5N^R^P^k(9+4N^) 

Also H.H.S. of Eq. (2.3.21) 

= rS ^( P ^^ + k ( 4 ^^( 9 + 4 N ^)"'’^^ 


- 1/2 


2Py\k O 

P^(-J^ - (9+4N^) ) 


Since 


dk 

31 ^ 


2k 

IT 


- 1/2 


R 


= 2 P ^ kR ^- 5 N ^ kP ^ R ^( 9 + 4 N ^) 


- 1/2 


. 3 . 21 ) 
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• (2,3,21) "therofore becoines, 

4R^P^+4R^P^k-5N^ll^(9+4N^) kR^=0 

k(2P^R^)= 5N^R^(9+4 n^) ^^^_4R^P^ 

—1 /2 

5n^(9+4n^) -4p, 

^ ....(2.3.22) 

Back-substitution of the values of k and P., into Eq. (2.3.22) 
gives 


* 5N^(9+4n^) ^^-4p., 

^/B,+6 = 5 


2P: 


-'I /2 1 /2 

5N^(9+4n^) -4(-3/2+1/2C9+4n^) ) 

2( -3/2+1 /2(9+4n^)^^^)^ 




1/2 1/2 

5N^+2(9+4N^) (3(9+4N^) -1 ) 

T72 T75 2" 

2(9+4N^) (1/2(9+4N^) -3/2) 
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APmffiI)C-D 


f 


From equation (2,4.20) on page number ( 27), we have, 

NR^ I.(4/3n) 

qa !— 

Iq (4/3N) + KNI (4/3N) 

•■applying the optimization criterion = 0, we get 

( (4/3i:) KNI^ (4/3N) ) ( nR^. (4/3N) ) 

= NR^I^ (4/3N).|^ (I^(4/3N)+ KNI^(4/3N),) 

; (2.4.21) 


L.H.3. of Eq. (2.4.21) 

= (I^(4/3N)+KNI^(4/3N))-(NR^ ( I^(4/3N) - (4/3N) . ) 

+ I^(4/3N) 4NR^+R^("|^ ) )) 


= (I^(4/3N)+KNI^(4/3N)).(NR^ ( ^ I^C4/3N)- ^ I^(4/3N) ) 

+1^ C4/3N) ( 4nR^-5/2NR^ ) ) 

= (Iq(4/3N) + KNI^(4/3N)) ( 4NR^I^ (4/3N) - N^R^Iq(4/3N) ) 

2 

= 4NiAq( 4/3N) I^(4/3N) + 4kNV (I^(4/3N) ) 

■_ ^ nV (Iq( 4/3N) )^* ^ kNVlQ(4/3N)I^C4/3N) 
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md R.H.d. or ITq. (?../+, 21 ) 


4 

n') “ / 


(I^(4/3 N)) ((_ (4/3 >:)+i^(4/3 

+ kN((I^(V3 N). N))C-m) 


NR‘*I^(4/3 n)(- ^i^( 4/3 t,)+ 5^ I^(4/3 N) 


5nk. 


lOKN^^ 


-^I,(4/3 N).i^I^(4/3 N)) 

nV(I^(4/3 N)) +2n2u\(I^(4/3 N)) - 

(4/3 N)I^(4/3 N) 


Lquating L.H,S* and i'^quation (2#4«21) can now be wrilrten 

as , 

4NI^(4/3N)I^(4/3 N) + (4/3 N)+4kN2(l^ (4/3N) )^- ^N^(I^(4/3 N))^ 

2 2 

+ •^N^(l^(4/3 N)) -2N^k(I^(4/3 N)) = 0 

2 

^N^((Io(4/3 N))2-(I^(4/3 N)) )-4NI^(4/3 N)I^(4/3 N) 
2N^(I^ (^i/3 N))^ 

substituting k = 1 /Bj^+6 ^ yields 

1 5 (Io(''/3 N))^ 5 2 

^“^(I,(4/3N))^ qwiHr-" 




FIG.I(a) LONGITUDINAL FIN-OF CONCAVE PARABOLIC PROFILE. 



FIG. Kb) CROSS SECTIONAL VIEV/ . 
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FIG. 2(a) LONGITUDINAL FIN OF CONVEX PAR>©0UC PROFILE. 



RG. 2(b) CROSS SECTIONAL VIEW. 
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FIG.3(a)SPlNE OF CONCAVE PARABOUC PROFILE. 



no. 3(b) CROSS SECTIONAL VIEW. 
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FIG. 4 (a) SPINE OF CONVEX PARABOLIC PROFILE. 



FIG. 4(b) CROSS SECTIONAL VIEW . 
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FIG. 5 OPTIMUM DIMENSIONS FOR CONCAVE PARABOLIC 
FIN. 


Fin Efficiency, -vv , Dimensionless . 



FIG. 6. FIN EFFICIENCY OF CONCAVE PARABOLIC FIN. 
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FIG. 7 OPTIMUM DIMENSIONS FOR CONVEX PARABOLIC FIN. 



F i n Ef f icie ncy, , Dimensi o nless . 
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FIG. 8 FIN EFFICIENCY OF CONVEX PARABOLIC FIN. 








Fin Efficiency, , Dimensionless . 







FIG.11 OPTIMUM DIMENSIONS FOR CONVEX SPINE. 


Fin Efficiency, , Dimensionless 





FIG. 13 OPTIMUM CURVES FOR VARIOUS FINS, 6=0. 



FIG.14 OPTIMUM CURVES EOR VARIOUS FINS, 8=1 
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36h 


2haL^ 


32h 


B] = 


28h 


% 

6 = 


S 20\- 


16r 



1. Concave Parabolic Spine. 

2. Convex Parabolic Fin. 

3. Concave Parabolic Fin. 

4. Convex Parabolic Spine. 








FIG. 16 OPTIMUUM CURVES FOR VARIOUS FINS, 6=3. 



Fin Efficiency . Dimensionless. 




Fin Efficiency , Dimensionless. 





Fin Effici 



V 

FIG.19 EFFICIENCY OF VARIOUS FINS, 'B| *6=2 



Fin Efficiency, >. ^ Dimensionless. 



Fie. 20 EFFICIENCY 
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1.4 1.6 


■q;? oTS bt? 0^8"” ^.0 

N — ► 

FIG. 21 COMPARISON OF OPTIMUM CURVES WITH OTHER PROFILE 

(fe'SO) 


Fin Efficiency, , Dimensionless. 



FIG.22 COMPARISON OF EFFICIENCY WITH OTHER PROFILES. 

{VBi*fe=0 ) 
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PROGRAM LISTING 


C 


i 

1 

1 

1 

1 

1 


13 


12 


22 


32 


42 


.S2 


PROGRAMME FDR OBTAINING THE OPTIMAL DIMENSIONS OF FINS 
REAL K,L,N 
integer FINNOM 

COMMON/ BE. SSEL/BFM(S00) ,BSM(S00) ,BFN<S00) ,BFF(£>O0) 

OPEN ( UN 1 1 '= 1 , F ILE= ' OPTD I M . OUT ' ) 

OPEN ( UN n =2 , F I LE= ' BFMBBM . INP ' ) 

0PEN(UNn=3,F ILE='BFNBFP. INP ' ) 

OPEN(UNIT=4,FILE='VALUES. INP') 

WRITE<*,3) 

FORMAT <//////////5X, 'PLEASE INPUT ONE OF THE F0LL0IN6 NUMBERS 
SPECIFY '/5X, 'THE TYPE OF FIN' // 

10X,'l =• CONCAVE PARABOLIC FIN'// 

10X,'2 = CONVEX PARABOLIC FIN'// 

10X,'3 « CONCAVE PARABOLIC SPINE'// 

10X,'4 = CONVEX PARABOLIC SPINE'/// 

40X,'THEN PRESS RETURN TO CONTINUE ') 

READ <*,*) FINNOM 
READ<2,*) <BFM(I> ,1=^1, 599) 

READ<2,*) <BSM<I) ,I«1 ,S99) 

READ<3,*> <BFN(I) ,I»l,fcO0> 

READ(3,*> <BFF(I) , 1=1 ,600) 

READ ( 4 , ) K , HF , HA , TF , T A , BELT A , QBI VEN 

PI«3. 1416 


L*» 0 . 0 
L**L +0-001 
Bs»K./ (HF*L) 

N»PNDN<:B, FINNOM) 

H*=-400*N 

W= (2 - «HA*L-»*2) / {K*N*«2) 

X=<K*TA) /L 
THETAF=<TF/TA-1. ) 

GO TO <12,22,32,42) ,FINNOM 
p«-0. 5+0. 5+SQRT ( 1 .+4*N^f*2) 

Q= < X*P*THETAF*W) / < 1 - +P* <B+DELTA/L) ) 
- ETA=’P/ < <N+-«2) * ( 1 .0+P* (B+DELTA/L> ) ) 


WRnE(*,*)Q 

( x2t^TAF.N.W*BFP > / (EFN («> i"> ‘ ’ 

eTA=BFP <M) / <N* <BFN (M) +N*BFP <M) * (B+DELTA/L) ) ) 

WRITE(*,*)0 
GO TO 52 

P=-J.5+0.5*SQRT<9.+4^tN**2) ,, 

Q:= < X*P*THETAF#PI*W**2) / < 1 - f 

ETA=-<P«3. >/< ^ 

WRITE<*,*)Q 

GO TO 52 /-vv / /ciCM/MV ^B+DELTA/D > 

0- ( X,lHETAF.ESH<m,N*Pl.W«*) / B 

ETA" < < 1 - 5/N) » <BSM (M) 3 ) / < <Br-n \ i 
WRITE<*,*>Q 


GO TO 52 
IF (Q.LT.QGIVEN) 
IF(FINN0M.LE.2) 
1F<F INN0M.6T-2) 
IF<FINN0M.LE.2) 


GO TO 13 

WRITE<1 ,5)L,W,N 
WRITE (1 ,7)L,W,N 
WRITE(1,9)Q,FTA 


ID 



1 


7 

9 

11 

C 


12 

10 

42 

52 

13 


62 


20 

72 

32 


30 

82 

C 


1 


1 

1 

1 


1 


IF(FINN0M.GT.2) WRlTEd , 1 1 IQ.ETA 

FORMAT (BX, 'LENGTH OF THE FIN = 'F7. 3, 'm ' 7/8X , ' W1D7 M [t n 
= ‘F7. 3, 'm' //ex, 'DIMENSIONLESS PAF^METER, N - ^ f; im 

FORMAT (8X, 'LENGTH OF THE SPINE « ■F7. 3, 'm ' //8X , 

SPINE «'F7. 3, ■m'//8X, 'DIMENSIONLESS PARAMETER, N 
FORMAT (/8X , 'HEAT DISSIPATION OF THE FIN = 'F10.4, ' ^ 

/8X,'FIN EFFICIENCY = 'F9.6) ^ 

FORMAT (/ex, 'HEAT DISSIPATION OF THE SPINE ='Fe.4' *^*RE'/ 

/8X, 'SPINE EFFICIENCY = 'F9.6) ' 

STOP 
END 

FUNCTION FOR CALCULATING THE VALUE OF 'N'- 
FUNCTION FNDN(B,FINNOM> 

INTEGER FINNOM 
REAL N1,N2 

CONMON/BESSEL/BFM(600> ,BSN(600) ,BFN(600) ,BFP(6»0O> 

GO TO (12,22,12,32) , FINNOM 
DO 10 1=1,600 

BINV=FNDBIN (0,0025*1 , FINNOM) 

IF (BINV.LE,B>GO TO 42 
CONTINUE 

IF (BINV.NE.BIGO TO 52 
FNDN=0. 0025*1 
RETURN 

N1 = .0025*(I-1. ) 

N2=. 0025*1 
FNDN=<N1+N2) /2.0 

IF (FNDBIN(FNDN,FlNNOM) ,GE.B) GO TO 62 
N2=FNDN 

IF ( (N2-N1>.GT- 1.0E-04> RETURN 
GO TO 13 
N1=FNDN 

IF ( (N2-N1 ) . GT. 1-0E-04>GO TO 13 
RETURN 

DO 20 1=1 ,600 

B3NV=<2, * (BFN (3 ) **2-BFP ( I )**2> - (BFP (3 > *EFN ( I > / < I* 
/(BFP(1)**2) 

IF (BINV.LE.B) BO TO 72 
CONTINUE 
FNDN=0, 0025*1 
RETURN 

DO 30 3=1 ,599 

BINV= ( (5. /3. >* (BFM( I >**2/BSrT( I ) **2-1 .0) ) - ( (2. *BFr>| 

/ (BSM(I)*(. 0025*1 +.0025) ) ) 

IF (BINV.LE.B) GO TO 82 
CONTINUE 
FNDN=0.0025*I 
RETURN 
END 

FUNCTION FDR CALCULATING 1/Bi 
FUNCTION FNDBIN (N, FINNOM) 

INTEGER FINNOM / 

REAL N 

GO TO (12,22,32,42) , FINNOM 
P=SQRT ( 1 . 0+4*N*»2 ) 


2 


42 


1 /(2*P*^i-2) 

RETURN 
RETURN 
END 



